0^ 



00 
O 
(N 

O 



X 



Variational Principle for Relativistic Fluid Dynamics 

Hans- Thomas Elze^, Yogiro Hama^, Takeshi Kodama^, Martin Makler'^ and Johann Rafelski'* 
^Instituto de Fisica, Universidade Federal do Rio de Janeiro, 

CP 68.528, 21945-970 Rio de Janeiro, Brasil 

^Instituto de Fisica, Universidade de Sao Paulo, CP 66318, 

05389-970 Sao Paulo, Brasil 

Centro Brasileiro de Pesquisas Fisicas, Rua Xavier Sigaud 150, 

22290-160 Rio de Janeiro Brasil 

Physics Department, University of Arizona, Tucson, AZ 85721 USA 

The variational principle for the special and general relativistic hydrodynamics are discussed in 
view of its application to obtain approximate solutions to these problems. We show that effective 
Lagrangians can be obtained for suitable ansatz for the dynamical variables such as density profile 
of the system. As an example, the relativistic version of spherical droplet motion (Rayleigh-Plesset 
equation) is derived from a simple Lagrangian. For the general relativistic case the most general 



^y^ • Lagrangian for spherically symmetric systems is given. 

-)— » ' 
CJ ■ I. INTRODUCTION 

O: 

I ' First applications of relativistic hydrodynamics to the process of multiparticle production in high-energy hadronic 

collisions can be found in the works of Fermi and Landau in the early 1950's 10,0 ■ Recently, extensive studies of the 

relativistic motion of fluids have been done with respect to the analysis of relativistic heavy-ion collision processes 

^ ' P~pi- III fact, a hydrodynamic description of high-energy hadronic and nuclear collisions has been successful in 

reproducing global features of these processes, such as multiplicity and transverse energy distributions. From the 

theoretical point of view, however, the foundation of the hydrodynamical picture for these processes is not a trivial 

matter. This is because the use of hydrodynamic equations of motion assumes implicitly the local thermal equilibrium 

via an equation of state of the matter. This means that the relaxation time scale and the mean free path should be 

^\ ,. much smaller compared to, respectively, the hydrodynamical time scale and spatial size of the system. In this sense, 

0^ ' one may wonder whether these conditions could easily be met for hadronic and nuclear collisions (for the collision of 

heavier nuclei they are expected to be approximately fulfilled for some specific scenario) . 
O ; On the other hand, from the kincmatical point of view, apart from the use of the equation of state, the equations 

I , of hydrodynamics are nothing but the conservation laws of energy and momentum, together with other conserved 
jr~^- quantities such as charge. In this sense, for any process where the dynamics of flow is an important factor, a 
(-H ' hydrodynamic framework should be a natural first step, at least at the level of phenomenology. The effects of finite 
• • I relaxation time and mean-free path might be implemented at a later stage by using an effective equation of state, 
incorporating viscosity and heat conductivity, or some simplified transport equations, see Ref. [pj and references 
therein. 
5J] ' Another important arena of extensive application of relativistic hydrodynamics is found in cosmology and high- 

er I energy astrophysics, such as the gravitational collapse of a stellar core to form a neutron star or a black hole, 
relativistic blast waves for the models of gamma ray bursts, etc. ||7|-p^. In these cases, the assumption of the local 
thermodynamical equilibrium is considered to be well justified. However, in the astrophysical applications we have not 
only to face the large scale systems but also to deal with the long range gravitational field simultaneously. For these 
reasons, the computer simulations of hydrodynamical scenarios for astrophysical problems usually become extremely 
expensive. 

The relativistic hydrodynamics is a local description of the conservation laws, written in terms of the energy- 
momentum tensor as 

a^T'^" = 0. (1) 

This is a set of coupled partial differential equations which, in general, arc difficult to solve exactly. Except for a 
few analytical solutions known for special cases, we have usually to resort to numerical solutions even for a simplest 
geometry, like one-dimensional or spherically symmetric cases. In the most of the cases, the numerical approach 
together with a realistic equation of state becomes prohibitively expensive, especially when coupled to some transport 
equations such as those for neutrinos in the case of stellar collapse or supernova explosions mW . 

In addition to the difficulties of solving the hydrodynamical equations, frequently we encounter with the situation 
where even the equation of state of the matter is not known precisely. Rather, we apply the hydrodynamical models 
to infer the properties of the matter involved in the process. In such cases we do not need the very precise local 



features of the hydrodynamical motion (for example, sound ripples, small local perturbations, etc.) but rather the 
global flow motion which characterizes the dynamics of the system assuming a given equation of state. 

For the reasons cited above, in spite of the presence of highly sophisticated techniques for hydrodynamic numerical 
calculations, some problems require rather simpler approaches which allow to analyze the dynamics of the system 
more effectively. In such cases, extremely local properties should be smeared out effectively, in order to extract global 
features of the flow more directly. As an example, we know that some global features of the high-energy hadronic and 
nuclear collisions can already be described by a simple fireball model p^ ]. Presently we aim at a dynamical scheme 
which improves the simplest fireball model in the direction of a more complete hydrodynamical description. 

We here introduce the method of effective action based on the variational principle to hydrodynamic equations 
of motion. As is well-known, the variational approach has practical advantages besides its formal side. Once the 
variational principle is established, we can use the method to obtain the optimal parameters of a given family of trial 
solutions. The effective Lagrangian and variational approach [^,|^, introduced to incorporate the effect of local 
turbulent motion in an effective way into a supernova explosion mechanism, is such an example. It was shown that 
such an approach is also useful to discuss the dynamics of a sonoluminescencing bubble in a fluid ||l^,|l^. There, the 
effective Lagrangian method was shown to be very useful in generalizing the so-called Rayleigh-Plesset equation to 
include in a simple way the effects of gas dynamics inside the bubble. There exist some analogous problems to the 
dynamics of a sonoluminescence bubble in the domain of relativistic energies, such as QGP or astrophysical fireballs. 
Thus, the relativistic generalization of the classical Rayleigh-Plesset equation will be useful. 

In the present work, we generalize the effective Lagrangian method to the relativistic hydrodynamics. By suitable 
parametrizations of the density profile of the system, approximate but very simple solutions of relativistic hydrodynam- 
ical models can be derived in this approach. In particular, we derive a relativistic generalization of the Rayleigh-Plesset 
equation and discuss the effect of relativity for the homologous motion of gas and fluid. 

Frequently the flow of the matter accompanies the production of entropy. In particular, when a shock wave is 
generated, the violent dynamical change of the density leads to a highly turbulent regime, which cascades into a 
smaller scale complex fluid motion and ultimately thermalize. In order to incorporate such effects of non-adiabatic 
pro cesses and simulate the dynamics of shock wave as a thin domain of non-adiabatic flow, Neumann and Richtmyer 
[ pTt introduced the method of the pseudo- viscosity which is still used extensively in many areas. We show that this 
approach can well be incorporated in our formalism and consequently, the relativistic generalization of the pseudo- 
viscosity method is easily obtained in our context. 

In the astrophysical applications of relativistic hydrodynamics, the inclusion of the gravitational field is essential. It 
has been discussed by several authors that the general relativistic hydrodynamical equations can also be derived from 
the action principle |18|-pi||. In this paper, we derive a simple general relativistic effective Lagrangian for spherically 
symmetric systems and deduce explicitly from it the equation of motion of Misner and Sharp |23] for gravitationally 
collapsing object. We also show that the concept of general coordinate system allows us to use a comoving Lagrangian 
frame in obtaining the effective Lagrangian of the special relativistic hydrodynamics. 

We organize this paper as follows. In Sec. II, we first review the variational formulation of the relativistic hydro- 
dynamics. Then, in Sec. Ill, we apply the variational scheme to spherically symmetric cases and establish an effective 
Lagrangian for the variational parameters of the density profile function. In the case of a homogeneous gas bubble 
in an infinite fluid, this equation is a relativistic generalization of the Rayleigh-Plesset equation well-known for fluid 
acoustic theories. In Sec. IV, we discuss nonadiabatic process and generalize the pseudo- viscosity of Neumann and 
Richtmeyer in the context of our relativistic variational principle. In Sec.V, we extend our approach to the general 
relativistic case, where, the metric functions are chosen as the dynamical variables. In the case of spherically symmet- 
ric system, we derive explicitly the effective Lagrangian for the comoving coordinate system. We discuss the relation 
between the comoving frame and space-fixed coordinate system. In Sec. VI, we summarize our present work. 

II. VARIATIONAL APPROACH 

Although not commonly found in general textbooks, the variational formulation of hydrodynamics has been studied 
by several authors |lq- |2l|j2J| -|26|| . For the sake of later discussion, let us first review briefly how the relativistic 
hydrodynamical equations of motion are derived from a variational principle. In the following, we take the velocity 
of light is unity, c = 1. Let the velocity field of the matter be 

v = v{r,t). (2) 

In order to keep the manifestly covariant notation, we express the flow in terms of a four- vector, w^(a;), where 

u° = 7, u = "fv. (3) 



and 

Uf,u^ = 1. (4) 

The flow of matter induces a cliange in the specific volume occupied by the matter. In order to facilitate the following 
discussion, we consider the case where there exists some conserved quantity, say the baryon number. Let the local 
density of this conserved quantity in the comoving frame be n. Then we have 

d^ {nu'') = 0. (5) 

We also define the specific volume V as 

V^l (6) 

Let us write the energy of the matter in this volume as 

E = eV, (7) 

where e is the energy density. The assumption of local equilibrium leads to the validity of the thermodynamical 
relations, such as 

where S is the entropy of the matter in the volume and P the pressure. In terms of the energy density this implies 

' d£\ e + P 



dn, g 



(9) 



Now take the action. 



Im^ j d^x[-e{n)+^{x)d^{nu^)+ icWK?/^-l)|. (10) 

and state the variational principle as 

5Im = 0, 

for arbitrary variations in u^,n,^, and Q. Then, as we will show in the following, Euler's equation for the relativis- 
tic fluid motion can be derived formally from the variational principlq^. Note that the last two terms in Eq.([lO|) 
represent the constraints among variables u^ and n. As we see in the next section, for the practical usage of this 
variational approach, it is convenient to choose the parametrization of u^ and n in such a way that the constraints 
are automatically satisfied so that the Lagrangian multipliers ^ and C, do not enter into the calculation. 

The variations in ^ and C lead immediately to the constraints, Eqs.fl) and (H). Applying an integration by parts 
to the second term in Eq. ([lO|) , the action can also be written as 

hi = j d^x |-£(n) - nu^'d^a^) + ^Cl^^) ("^^m - 1)| ■ (11) 

We remark that the derivation in this section is equally valid in general coordinate systems. In this case, the partial 
derivative, d'^, in Eq.(ni|) should be replaced by the appropriate covariant derivative, and correspondingly in the 
following equations. Furthermore, the volume element d'^x should be replaced by the invariant volume element 
V^d'^x [|3. 

The variation with respect to n leads to 



^In fact, the action above only applies to the case of non-rotational flow. It is also possible to formulate the variational 
principle for more general flow pattern. See |19|,^ for discussion. 



5n 



— - M^a^e = 0. (12) 



Note that, if the motion of the fluid is not adiabatic, then Se/6n does not necessarily be equal to the usual derivative 
de/dn (see the later discussion). On the other hand, the variation in u^ leads to 



and substituting this into Eq.(|l2|), we obtain: 



^mC = ^C"m: (13) 



C = -"|^, (14) 

on 



where we have used u^ — 1, cf. Eq.(Q). Thus, Eq. ([l3|) is rewritten as 



df,i = -—Ufj,. (15) 



Taking the contraction of both sides with the four- velocity, we have 



.'^a,C = -|u,.^ = -|. (16) 

Thus we conclude that 

^ = -- (17) 

since u'^df^S, — d^/dr. 

In the usual hydrodynamic equations, we assume that the matter is always in thermodynamical equilibrium. Fur- 
thermore, if there is no viscosity or heat conduction, the energy change associated to the motion is adiabatic, that is, 
the change in specific energy E caused by the change in the specific volume V is given by 

SE = -P6V, (18) 

where V = 1/n and P is the pressure. In such cases, we have 

ds 5{nE) _e + P 
dn 6n n 

and 



(19) 



^Sn J n 
Therefore, for adiabatic changes of the density n, Eq.(|l7|) becomes 



d[^)-^ '-dP- (20) 



d^ 1 

On the other hand, we have also 



^MTT = — ^M^- (21) 

OT n 



e + P 






= ^^M^ - ^-^^M T^^uu'' = — 9,,C, (22) 



that is, the two derivatives 9^ and d/dr commute when applied to ^. Therefore, from Eqs.(|lq) and (17) we can 
eliminate ^ to get 



vJ'd, 



d^P 



Again using Eq.(EQ), this reduces to 

u^ u-{d,P) + {e + P) u-d, u^ = d,,P 
The first term of this equation is further modified as 

'e + P] 



(23) 



(24) 



u''{d^P) = u''nd^ 



= u'^d, (e + P)- 



P 



u'^d^n 



= u^d, (e + P) + (e + P) d.u'', 
where the continuity equation (^ was used. Thus, Eq.(|24|) can be rewritten as 



d.T-;^o, 



where 



T^i^ = (e + F) Uf_,u^ - Pgf_,„ 



(25) 



(26) 



(27) 



That is, we arrive at the equation of motion of relativistic fluid dynamics with the energ y-m omentum tensor of the 
perfect fluid. From this equation, we obtain the relativistic version of the Euler equation |271|, 



— u + [v -V) V 
at 



{e + P)j^ 



VP 



^dP 
'~dt 



(28) 



In the above derivation, we assumed the Minkowski space-time metric, but as mentioned before all the calculations 
can easily be extended to the case where the metric is more general. For example, the final result Eq.(p6|) in the 
curved metric is 






(29) 



where ; stands for the covariant derivative, as usual. 

As pointed out by several authors |19j2q,pq], the above scheme leads only to non-rotational flow. This can be seen 
from Eq. ( [l3| ) , where the velocity field is proportional to the four-gradient of a scalar function. In order to include the 
rotational flow, we have to add a term coming from another constraint with respect to entropy in the original action. 
However, for the spherically symmetric case below, one does not need to worry about the rotational flow, so we omit 
the discussion for the sake of simplicity. 



III. SPHERICALLY SYMMETRIC CASE 



The above variational approach is particularly useful when we can solve the constraint equations explicitly. One 
dimensional, or spherically symmetric system is such a case. Here we study the spherically symmetric case. Let the 
density profile in a space flxed frame be 



n = /(r,t). 
Then the velocity field is determined from the continuity equation as 

1 



r^f Jo 



(30) 



(31) 



If we use this expression for the velocity field, then we can omit the constraint terms in the action. Thus, we have a 
Lagrangian for / as 



L = -47r / r dr e (n) , 
Jo 



(32) 



where the local density in the comoving frame is given by 



7' 



(33) 



To see explicitly how the variational principle works with this Lagrangian, let us consider an arbitrary density variation, 

f^f + Sf. 
Under such variation, we get 



and. 



SI = -4tt I dt I r^dr 6e (-/ 



^Att I dt I r^dr (^\ (-^ + infSv ] , 



\5n J \ 7 



Sf 1/2, 

/ r^f Jo 



Sv = — —V 5— r / r 5fdr 



(34) 



51 ^ All / dt 







on 



'Ait I dt 
-Att / dt 



r'^dr i — ] ^5f 

^'"^{(Ir dt 



l-v-fi—Sfv 7T / r Sfdr 

dr —— ] vj I r' 6fdr' 



dr { — ] v^ 



51 



From this we have 



5e_ 



dr { -— ] V"f 
on 



= 0. 



Se\ 1 


d 


'f5e\ ■ 


on J 


dr 


\5n) _ 



Taking the derivative with respect to r of the both sides, we get 

d_ 
di 

Using the adiabatic relation of the energy density and pressure in the equation of motion, we get 

dv _ 1 (dP dP' 

which is again the relativistic Euler equation for spherically symmetric case 



(35) 



(36) 



A. Effective Lagrangian 

Now let us introduce a parametric ansatz for the density profile as 

n(t,r)=/(r,a(t)), 



(37) 



where a — a(t) is a (set of) time-dependent parameter(s) which determines the form of /. We suppose that / is 
suitably normalized, 

/•CXD 

Att / r^fdr = N. (38) 

Jo 



The velocity profile, Eq.(|31|) becomes 



-ir,t)=-a-^^ 



with 



oa , 



r'^('J-)dr'. 



A-A(r,a) = - — 
For more than one parameter, a = {a\i = 1, ...,n} Eqs.(pl|) and (Q) should be understood 



(39) 



(40) 



(41) 



The important point here is that the velocity field is linear in d. The effective Lagrangian for our variable a = a(t) 
becomes 

/•CO 

L(a, d)=-47r/ r^dre{n), (42) 

Jo 



where n = f{r,a)/j and 7 = ^/y^ ^ (^ ^) ■ The equation of motion for the variable a is obtained from the 
Euler-Lagrange equation 



d fdL\ dL 



dt \d: 



a I da 



= 0, 



which is written as 



For adiabatic motion, we get 



An j r dr 
/o 



d ( de 



dt 1 Sn 



n-y^A'^d 



Se dn 
5n da 



d 
dt 



/•OO nOO 

/ r'^dr{{e + P)j'^vA} = - r^dr{e + P) 

Jo Jo 



= 0, 



19/ o dv 

7 V 



fda 



da 



(43) 



(44) 



(45) 



The effective Hamiltonian H is then 

oa 



Se .9/1 
T-na-^ {- \ +e 
on oa \7 



= Att r dr 

= AiT I r^dr [{e + P)-f^ -P]= d^r T 
which is in fact the total energy of the system and a conserved quantity. 



47r / r^dr [{e + P) -f^v^ + e] 



j3„ rj^OO 



(46) 



B. Relativistic Rayleigh-Plesset Equation 

For an example, let us consider a system composed of a homogeneous spherical gas bubble surrounded by a 
homogeneous fluid. We then introduce the ansatz, 



O JSJ 



(47) 



where the radius of the gas bubble R = R{i) is the only dynamic variable. We use the subscript G and L to specify 
the quantities in the gas and the fluid, respectively. For example, Nq and Nl are number of particles (constant) in 



the gas and fluid, respectively. The outer (constant) radius of the fluid Roo is introduced here to take into account 
the conservation of the number of particles in the fluid, but ultimately should be taken equal to oo. The velocity fleld 
is then determined as 



v{r) = ^R, 



R^ (Rl - r3) 
r2 (Rl - i?3) 



0<r <R, 



R, R<r < Ro 



(48) 



so that 



AM ^ ^. 



0<r <R, 



R^ {Rl - r3) 
r2 (i?g^ - i?3) 



R<r < Rr, 



(49) 



The effective Lagrangian for R is then 



L^-^T^j r^dreG{.fG\jl-[^R] | - 47r 



ir-'^'A^-im^ 



R 



(50) 



The equation of motion for R = R{t) is given by 



d_ fdL\ _ dL_ 
di[dR) ~dR' 



(51) 



where 



^=47r— / drr'*(eG + PG)7^+47ri?^^ , 

dR R Jo Jr ''~ 



dr , _ . T / i?iL - r^ ^ 



ieL + PL)r 



2 / ^i-oo 



RL-R^ 



(52) 



and 



dR 



= -47ri?2 [ec - eL\R - 47r / r'^dr [sg + Pi 

Jo 



'-^|-|+^' I li 



flo 



-4./^ .^rfr(s. + P.) 3^_^3 



At this stage, we can take the limit R^ -^ oo. Thus, we have 



3i?2 ^R^2Rl + R^){Rl- 



{Rlo'R'f 



-R 



b2 



(53) 



aL 
aR 



7? f" f°° dr 

47r^ / dr r^ {eg + Pg) t' + 4^P^i? / ^ (^l + P. 
^ Jo Jr ''" 



'L)7^ 



(54) 



and 



— = -47rp2 ([sG - eLJR + [SL + Pl]^) + ^ J ^^dr {eg + Pg) 

-A^ prHr (eG + PG)7'[|] "^ + SttR^R^ J^ ^ {sl + Pl)i^. 



The equation of motion then takes the form, 



(55) 



d_ 
di 



{h + h)R*R =F-{h-2h)R'R\ 



(56) 



where 

h= f dxx^ieG + PG)l^ (57) 

/•OO 7 

l2= ^(sl+Pl)!'. (58) 

Jl X 

with a; = r/R and 

F = -i?2 ([£G - el];, + [SL + Pl]„) + j^J r^dr {ea + Pg) , (59) 

where the subscript represents the quantity evaluated at ti = 0. Eq.(pq) is a full relativistic equation of motion for 
the radius of a gas bubble under the homologous motion of the system. 

To see the non- relativistic limit of Eq.(pq), we separate the energy density into the sum of the rest- mass energy 
density p and the internal energy density e™* as 

e^p + e'"-*, (60) 

and expand Eqs.( |57|j58|j59| ) in a power series of small parameters such as , w^, ei„f/p,and P/p in the non-relativistic 
regime. We have 



Thus, 



and 



{e + P)-f^^Po + el"*+Po + ^Pov^ (61) 



Jl c. i {po.G + e^l + Po.g) + ^Po,gR^ (62) 

l2^PO,L+eo% + Po,L + ^Po,LR^: (63) 



F:^R^ (Po,G - Po,l) - Q {po,L + Po.l) + ^ {po,G + Po,G + c^sPo,g)] R^R^ 

-\ (pox + yPG,o) R^R^- (64) 



In the lowest order in R we get 

1 \ „;:. 3 



tPo,g + Po,L 1 RR + ^Po,lR^ = Po,G - ^0,L, (65) 

which is the usual Rayleigh-Plesset equation of a gas bubble inside a liquid ||l^ without the term for the energy 
dissipation due to the sound radiation. 

When the equation of state of the fluid and the gas are given as 

P(x p^ (66) 

where F is the adiabatic index, all the integrals in Eq.(^6|) are expressed analytically. We get 

MR^ii + 1^po,gJi + Po,lJ3 + ^^1^ (^" + P)o.G ^2 + (e" + P)o,l J^] ^'^' 
= -i?2 ([eg - eL]n + [sL + Pl]o) + R^ (po,g Jg + (e" + P)o^g "^lo) ' (67) 

with 



M = ipo,G Ji + P0XJ3 + \ (e™ + P)o G ^2 + (e" + P)^j^ Ji 



\ {po,gJ5 + (2 - r) (e- + P)„_g Je) + i (po,L J7 + (2 - T) (e™ + P)^,^^ Jg) 



i?^ 



(68) 



and 



J3 = F 
J5 = F 
Jr = F 
J9 = F 



'1 


5' 


2 


2 


'1 


1' 


2 


4_ 



2' 4 'T'^'' h ^^4 = -F 

3 7 
2' 2 
3 5 
2'4 
3 1 
2'~2 



-,P" , J6 = F 



-,R ] , Js — F 



- , P^ I , Jin — F [ - , — — 



1 


r 5" 

~ 2' 2 


I..'). 


2'4 


•!«=)' 


2' 2 


4,«=), 


2'4 


■l^)^ 


V 


■3 r' 

2'~2 


.|«^), 



where F{[a, b], c, z) is the hypergeometric function. 

For the sake of illustration, we show in Figs.l, 2, and 3, time dependences of the radius and velocity described by the 
relativistic Rayleigh-Plesset equation Eq.(67) for 3 different initial conditions. In this example, we consider the case 
where the both gas and fluid have the same mass density and the adiabatic index F = 4/3. Three cases shown here 
are for the different values of the initial gas pressure, {P/p)^^ p = 1/10 (Fig. 1), 5 (Fig. 2), and 100 (Fig. 3), keeping 
the ratio of the initial gas to the liquid pressure Pg,o/Pl.o = 100. Solid lines are for the full relativistic equation of 
motion (Eq.O) and the dashed ones are for the non-relativistic limit, Eq.(p5|). For the low initial gas pressure, two 
solutions coincide (Fig. 1). For the extremely high initial pressure, the motion of the bubble becomes completely 
relativistic (Fig. 3) and the non-relativistic equation of motion differs completely from the relativistic equation. Note 
that in the relativistic equation, the velocity saturates at v/c = 1. Of course, in this extreme example, the fluid motion 

becomes supersonic (v>Vs — yl/S = 0.577] and the hypothesis of homologous motion may breakdown. However, it 

is important to note that, there exists a case where the non-relativistic approximation fails down completely although 
the fluid and gas motion are still subsonic like in the case of Fig. 2. 



IV. NONADIABATIC PROCESSES 



In many cases, the change of the density associated with the flow of the matter causes non-quasi static processes 
within the hydrodynamical volume element established in the practical calculations. For example, in the limit of large 
Reynolds number, the dynamical change of the volume easily leads to a highly turbulent regime in small regions of the 
fluid, and this complex fluid motion will gradually thermalize inside the volume element. In such a case, there appears 
the heat production inside of such a volume element. If the time scale for the thermalization is negligible, then the 
heat production can be expressed in terms of viscous tensor, and the hydrodynamic equation of motion becomes, 

df, (T^'' + E^") = 0, (69) 

where S^'' is the shear tensor. When there is no heat transfer, we may take 

1:'"' =q{u''u''-g'"'), (70) 

where q is a function of the local thermodynamical quantities, like p, P and its derivatives. From Eq. (|69|) we have 



U^d^,(^^)+iP + q)u^dJ^] =0. 



(71) 
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This means that, in the Lagrange comoving system, the change of the specific energy with respect to the proper time 
is given by 

dE dV dV 

Thus we identify the function q as the rate of the production of the entropy S with respect to the volume change, 

rr^dS dV 

The specific form of S^'', Eq.(^Q) allows us to write 

f^"" = {e + P + q) u^'u" -{P + q) <?^^ (74) 

which conserves 

From the above conservation law we get immediately the equation of motion, 

dv 1 (d{P + q) d{P + q)\ 

which describes the relativistic hydrodynamical motion under the local entropy production, Eq.d72). The function 
q should be specified appropriately according to the non-adiabatic processes representing the conversion of kinetic 
energy of the collective motion to the internal energy of the matter. Such a viscosity was first introduced in the non- 
relativistic hydrodynamics by Neumann and Richtmyer ||17| in order to simulate the entropy production mechanism 



at the shock front. Eq.(76) is the relativistic extension of the method of pseudo- viscosity of Neumann and Richtmyer. 



The above scheme is easily incorporated in the variational formalism. In the presence of non-adiabatic processes, 



the variation in the specific energy in the previous section, Eq.(18), should be replaced by 



SE = -P5V + (5Q = -PSV + T5S ^-{P + q)SV (77) 

where SQ is the generated heat associated with the non quasi-static density variation. Consequently we should, 
instead of the adiabatic relations (ll^JSO) , use 

de e + P + q 
on n 

'^(t) -^-d{P + q)- (79) 

\on J n 

From these substitutions, we get immediately Eq.(|7^). In terms of parametric representation, the equation of motion 
is given by 



±j\^dr{{e + P + q)^^v^]=-j\HT{e + P + q)(^-^J--^\^£ 



(80) 



The effective Hamiltonian for dynamical variable a is given again by 



oa 



Mr [{e + P + q) 7^ - (P + q)] = f d^r {T"° + Y.°") , (81) 



f.o, 



= in 
which is conserved. 



for the equation of motion, Eq.(BOp, together with Eq.(f72). This H can again be identified as the total energy of the 
system including the internal heat energy generated in the fluid. 
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V. GENERAL RELATIVISTIC HYDRODYNAMICS 

For the application of the present formahsm to astrophysical problems it is essential to include the effect of gravity 
through the theory of General Relativity. The variational approach of the general relativistic hydrodynamics has been 
discussed by several authors p"8|-pl[. In this section, starting from the variational approach, we show that the method 
of effective Lagrangian can also be established taking the metric as one of the variational trial functions. Let us first 
review how the general relativistic energy and momentum tensor are derived from the variational approach. 

A. Energy-Momentum Tensor and Einstein's Equation 

The total action is given as 

I — Ig + Im, 
where 

iG^^Jd^xV^U, (83) 

is the action for the gravitational field and k — SnG with G the gravitational constant. As usual, g = det Igui^ii is the 
determinant of the metric tensor (7^,^, and SR is the curvature scalar. The action of matter is now given by |18| 

.4^1 1 



hi = J d''xV^i-ein)+ax) (.nu^,^ + ^ax) {u^'u^ - 1)| , (84) 

where ^ and C,/2 are Lagrange multipliers as before. As usual, ";" represents the covariant derivative and the factor 
-^—g is inserted to guarantee that the Lagrangian density is a scalar. The variation of the action should be carried 
out with respect to (7^j/,n, u'^,^,and C, independently. The results of variations with respect to n, u^,^, and C are 
the same as before (see Eqs.(O), (0^,(0), (Iq) , (p9) and the comments for the covariant derivative in Sec. II). Thus, 
these variations gives the relativistic hydrodynamic equation for a given metric g^''' . 
The functional derivative with respect to g^i, is calculated to be 

5Im _ ^N/^{£(„)+„(a^^(^))yM}+ 1^^3^j,M^-^ (85) 



5gt,p dg^^ 2 

where we already employed the constraints, cf. Eqs.(^,||). Using 

'^- iv^.-, (86) 



and substituting the values of 9^^ and C, we get, 

^^'' = \^r~9r W) - n i'-^:^u^ u^] + \v~g{e + P) u^u^ 



Sg^i, 2 [ \ n J ] 2 

^ ^{{e + P)u^u" -Pg^"}. (87) 



2 
Comparing this result to the definition of the energy-momentum tensor 





Shi _ 

Sgtiy 


1 
'2^ 


/^T'^^ 


rjifiU 


= {e + 


P)' 


i^M" - F 



we identify that 

Pg^". (89) 

which is nothing but the energy-momentum tensor of the perfect fluid. Thus, the energy-momentum tensor of the 
fluid is derived from the Lagrangian density Eq.(p3) just as in the case of the fleld theoretical Lagrangian. Note that 
the role of constraints are essential for this derivation. 
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The variation of the gravitational action Iq with respect to g^^ gives the usual Einstein tensor Q^^ , 

SIg _ 1 



-gG^", 



so that we get 

gt"" ^ KT^"', (90) 

which is the Einstein equation, as expected. The hydrodynamical equation, Eq.(|2£ 

T^;; = 0, (91) 

can be re-obtained from Eq.(po[) due to the Bianchi identity, 

^'^ ;m = 0. (92) 

It is interesting to note that if we use the metric functions as basic dynamical variables then, the hydrodynamic 
equation of motion is obtained somewhat indirectly from the properties of metric tensor and the variational principle 
does not lead directly to the equation of motion. This point will be discussed later again in the context of the 
derivation of the special relativistic equation of motion using the comoving coordinate system. 

B. Spherically Symmetric System 

The derivation of the equation of motion above is too formal and not much useful to be applied directly for some 
practical problems. To make use of the variational approach, it is necessary to establish appropriate trial functions in 
order to write down the effective Lagrangian for these functions. As in the case of special relativity, this is possible 
when the system has appropriate symmetry, such as spherically symmetric distribution of matter. Many problems 
of the gravitational collapse of stars, the structure of neutron stars, and the Robertson- Walker cosmology can be 
discussed in this symmetry. Here we establish the effective Lagrangian for the spherically symmetric system. 

The most general form of the metric for a spherically symmetric system can be taken as ^7| 

ds^ = e^^dT^ - e^^de - r^dn^, (93) 

where (T,^) denotes the time and radial coordinates and (j) — (f){S,,T), A — A(^,T), and r — r{^,T) are unknown 
functions to be determined. Usually, if we consider the radial velocity field of the fiuid as an independent variable, 
then we need only two independent functions in the metric and we may choose, for example, r = ^. However, with 
the above metric involving three functions, we can further take the so-called comoving frame in such a way that the 
space-like components of the four- velocity field vanish everywhere [ p8| , 

m'^ = (u", 0,0,0). (94) 

From the normalization condition w^w^ — 1, we get 

u" = e-^. (95) 

In this comoving frame, the conservation law is expressed as 

(nu^).^^ = -La^ {V^nu^) = -^^Ot (r'e^n) = 0, (96) 



so that the density n of the conserved quantity, say the baryon number, is given by 

P 



(97) 



where p — p(^) should be determined by the initial condition. In this choice of the metric, the matter Lagrangian 
density is expressed as 

Cm = -A/^e(n) = -e'l'e^r'^e{n), (98) 
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where n is given by Eq.(g^). No terms with Lagrangian multipUers appear, because the constraints are automaticaUy 
satisfied. The gravitational part is calculated as 



Ca = 



1 

2k' 



e e^ 



r {2r(j) + r ) e 



i\ „-2A 



[2r\ + ; 



-20 



1 



(99) 



where we introduced the notation / = df /dT and /' = df/dS,. In the above, we omitted the part which can be 
written as the total derivative of a function, since this does not alter the equation of motion. The total Lagrangian 
of the spherically symmetric system is then given explicitly as 



C [(f), A, r] = e 



- ^'t'^^^2 



e r 



-e(n) 



1 



2Kr^ 



r' [2r(t)' + r') e 



'\ „-2A 



{2r\ + f\ 



-24> j^l 



(100) 



When the variation for functions 0, A, and r are in fact arbitrary, this Lagrangian is equivalent to Einstein's equations. 
To see this, we write the Euler-Lagrange equations of motion for 0, A, and r from this Lagrangian to get 



_ 1 



1 ,, fr' , r" r'2 



r 



2e-^ ( ^ + -A 



(101) 



and 



P 



-2e 



-2A 



J 2 



2r2 



2e 



-20 



2r2 



P 



2 A I /// I 1 12 il \l 



A' -f 1 (r" + 0'r' - A'r') 



(102) 



(103) 



;2 



W^ A -F A^ - 0A -F - f -h f A - (/)f 



We verify directly that these three equations are exactly those corresponding to the diagonal part of Einstein's 
equation. In fact, writing these equations in the form 



1 






0' 



(104) 



'-y 1-, 

K 



(105) 



and 



P 



-Q\. 



(106) 



wher e Q\ Q\ and Q\''s are defined, respectively, as the quantities in the square bracket [ ] of Eqs.(lOl), (102), and 
(103), we can identify the functions Q as the diagonal components of the Einstein tensor corresponding to the metric 

(Ef- . 

The only difference between our formalism here and Einstein's equation is that in the former there is no equation 
corresponding to the non-diagonal element, Q\ in the latter. In Einstein's theory, this quantity should be zero. 



^0 



2e 



-2A 



U' - r<j)' 



Xr' 



(107) 



This is because in the comoving frame the energy-momentum tens or T^' is diagonal. There fore , to pro ve that our 
result is identical to the usual theory, we have to show that Eq.(107) is a consequence of Eqs.(104) - ( 106 ). Although 
this proof is rather basic matter and could be found in text books, we show it explicitly for the sake of later discussion. 
We first start with the well-known Bianchi identity (for example, see [g^, p. 363), 



Gt ;m - d, {V^Gt) + \^^[gc. pd^g'^^] = 0. 



(108) 
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In our case, the first component f — leads to 



do {V^G"o) + di iV^Gl) + ttV^ 



E.9-^"a5o5"" 



On the other hand, we have 



do {v^g°o) 



a=0 



^gaaG^adog' 



.a=0 



r 



r^e^e- ^ ^ + 2^ {g\ - g\) + {G^ + g\) [2'- + \ 



r^e^e" 



dt 



£ + (£+p) 2- + A 



where Eqs.(104) - (106) together with the energy conservation. 



de d {pi 



r^e^l 



dn dT 



--(e + F) 2- + A , 



are used. Comparing Eqs.(109) and (110), we get 

di iV^g'o) = 0, 



(109) 



(110) 



(111) 



2^0^A^l _ 



e^e^^l, = CiT). 



(112) 



where C is a function of T only. For a non-singular metric we should have r(^ — 0,T) — 0, hence we conclude that 
C(T) = 0. Therefore, we obtain 



(113) 



g'o - 0. 

This completes the proof that our result is equivalent to Einstein's equation. That is, the Lagrangian density Eq.(|lO 
describes correctly the dynamics of a spherically symmetric system of an ideal fluid and gravitational field 



C. Misner-Sharp Equation 



Together with Eq.(107), Eq.(113) implies the following relation, 

f' - f(p' - Xr' ^ 0, 



(114) 



which can be obtained from Eqs.(104) - ( |106| ) directly, without referring to the Einstein tensor G^ and its properties. 
Following Ref. 123] we can express the equations of motion in a more convenient form. First, putting G^i = in 
the second component of the Bianchi identity, we have 



dg\ 

dr 



+ 2^-{g\~g\)^-{g^, + g\)<t>'. 



Substituting Eqs.( |l04| ) - ( |106| ), we get immediately that 

P' = -(£ + P)0', 

which is the Euler equation in the comoving frame. Now we introduce a quantity U defined by 



dr 



(115) 



(116) 



(117) 



where r is the local proper time and d/dr is the total derivative. The relation (114) is expressed in terms of U as 
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^0; U' 
e ^X — —-. 



Now Eq.(lOl) becomes 

STrGer^ = I - rr' [e'^^)' - e'^^ (r" R + r'^) + U^ + 2rUU'r'~\ 
which can be integrated as 



-2A 



rf2 



l + C/2 



2MG\ 



where, as before, G is the gravitational constant and 

/o 



M{^,T) = 47r / er^r'd^ = Att er^dr. 
Jo Jo 



Inserting Eq.(119) into (105), together with (116), and after some manipulations we obtain: 



dr^ 



Ue 



-0 



e + P 



1 + U^ 



2MG\ fdP 



dr 



M + AttWP 



(118) 



(119) 



(120) 



(121) 



This form of the equation of motion was first obtained by Misner and Sharp ||2^. Equations (116), (117), (120) and 
(121) together with an equation of state completely determine the dynamics of a spherical collapse, or bounce, which 
might be relevant for the study of the gamma ray bursts Jl^ . 



T o see the relation between this expression and the special relativistic Eulcr equation Eq.(|36|), we set G = in 
Eq-dH), 



d'^r 
d^ 



(^)(--^)(f). 



Let (t,r) the coordinate system fixed in the Minkowskian space-time. Thus, the line element is given as 



ds^ = dt^ - dr^ - r^dfl^. 



in a space-fixed coordinate system. However, in the comoving coordinate system, we need non-trivial metric functions 
as follows. First we introduce a coordinate transformation. 



where we may identify the coordinate ^ as the comoving Lagrangian coordinate. Thus, by definition. 



dt 



(122) 



(123) 



^^cont. 



It is always possible to choose the variable T so that the cross term in the above equation vanishes and the line 
element can be written in the form of Eq.(|93|), 



Here, 



ds^ = e^-^dT^ - e^^dC^ - r^{£,, t)dn'^. 



dr = e^dTU^=o, 



(124) 



(125) 



is the (local) proper time. Note that a local Lorentz transformation relates the infinitesimal coordinate differences to 
those of proper time dr and local radial distance e'^d£_ by 



dt 
dr 



7 vj 



dr 
e^d£, 



(126) 
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On the other hand, since 



we conclude 



In this way we have 






(127) 



(128) 



dr\ d f dr 

dry dfj^^^'diy'dtj^ 






^7 



9t;^ 



Now, from Eq.(126) we have 



7Ml + 7V)(|)^ 



dr\dT=o = le^d^, 



(129) 



so that 



dP 
dR 



dP 
dP" 



dp\ fm 

dt)Xdr 
dP 
'dt 



We also have 






5^y/ w* Wc " vs^y/'' w/ 



(130) 



(131) 



Therefore, Eq.( [l2l|) becomes 



f dv\ f dv 



1 1 



^ e + PY 

which is exactly the relativistic Euler equation (|3q). 



dP_ 
dr 



dP 

'at 



D. Variational Principle in Comoving Coordinate for No Gravity Limit (G -^ 0) 

The above discussion suggests the possible use of the comoving (Lagrangian) coordinate system even for cases with 
no gravity, that is G ^ 0. The effective Lagrangian presented in Sec. Ill is based on an ansatz for the solution of the 
continuity equation in the space-fixed coordinate system. By using a comoving Lagrange coordinate system, we may 
better choose the trial function on physical grounds. Of course, the two systems of coordinates, in principle, should 
be equivalent if the ansatz has enough flexibility to the express any arbitrary flow pattern of the matter. However, for 
practical applications, the appropriate choice of the coordinate system is essential to get better results. For example, 
it is technically difhcult to introduce the velocity dependence in the ansatz for the density profile consistent with 
the continuity equation. Therefore, for an ansatz like Eq.(B^, established in the space-fixed coordinate system, the 
relativistic kinematical effects may induce some spurious effects on the dynamics of the parameters. On the other 
hand, if we can choose the parametrization in the comoving coordinate system, such kinematical effects are expected 
to be automatically included in the equation of motion. 

In the limit of G -^ , the space-time reduces to that of Minkowski and obviously the gravitational part of the 
Lagrangian density (100) vanishes. However, in the comoving frame the line element has still the form (p3|) and 
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functions e'*' and e remain unknown. If we drop out the gravitational part from the action, the variati onal principle 



does not give information on these functions. What should be done in this limi t is that Eqs.(116) and ( |119| ) are used 



as constraints among the unknown functions, r, (j) and A. Setting G = in Eq.(119) and using the relation 



we obtain 



e^ = -. (132) 
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On the other hand, assuming the isentropic initial condition and adiabatic process, we can integrate Eq.(ll(;) with 
respect to ^ to get 

e-^ = ^^ = h{n) (133) 

mn 

where h is the specific enthalpy of the matter and m is the rest mass of the constituent particles. This integration 
constant was chosen so that e"^ ^- 1, n ^- 0. 
Now the action becomes 

I = - I dT I d^e^e\''e{n) = - f dT f di—-^e{n), (134) 



where 

and p = p{^) is determined from the initial condition. 
Let the T-dependence of r be specified as 



(135) 



r{^,t)=f{^,a{T)). (136) 

Then the Lorentz factor 7 is expressed as 



2 
d2 



7 = Vl + e-2^r2 = I + /j2(„) f^\ 
so that the number density n — n{a, a; f ) should be determined by the equation 



f'f'n^pi0Jl + h^n)(^^^ a\ (137) 

Finally the effective Lagrangian for a — a{T) is given by 

L{a,a) = - j dip{^)e{n). (138) 

The Euler-Lagrange equation of motion for a then takes the form 



(139) 



Differently from the case of the space-fixed coordinate syst em, we need the equation of state to determine the 



density of the matter as a function of our dynamical variable, (136). Note the difference between the density profiles 



ri(r, t) of the Sec. Ill and n(^, T) defined by Eq.(137). The former is defined for constant time, t = const, of the space 
fixed global coordinate system, and the latter is defined for constant time coordinate, T = const, of the comoving 
coordinate system. 
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Although the comoving coordinate system seems physically more advantage ous than the space fixed global coordi 



nate system in choosing an ansatz, it may generate a difficulty in solving Eg. ( 137 ) for a given equation of state. For 
the ideal gas, like 

P ex n'*' 



with 7 = 5/3 or 4/3, an analytic solution of Eq.(137) for n can be obtained explicitly. However, for general cases. 



analytic solution is not available. To avoid this, one may be tempted to use the proper time defined by 

dr = e'^dT (140) 

instead of the time coordinate T and introduce the ansatz, 

KC,r)=/(^,a(T)), (141) 



in substitution for Eq.(136). In this case, the Lorentz factor becomes 



->/'+ii{)'(sy '-' 



involving no density dependent term in it. Thus the density is expressed directly as 



fhhmm 



without need for solving the equation, (137). However, unfortunately, the pair of variables, (t, ^) do not constitute 



the proper integrable coordinate system, so that the boundary condition for the variation principle on the action, 

I = - I dr I d£.e^r^e{n) (144) 

is not properly defined and a simple Euler-Lagrange equation for fixed r leads to a wrong result. 

VI. DISCUSSION AND CONCLUDING REMARKS 

The variational approach for systems of fields, including the general theory of relativity, is of course a basic and 
standard theoretical framework and has been well studied, even for the application to the hydrodynamics discussed 
here. However, to the authors' knowledge, except for the formal derivation, no explicit variational formulation for the 
practical application for relativistic hydrodynamical systeins has ever been carried out. 

From the formal point of view, matter described by hydrodynamics is rather a phenomenological concept than the 
consideration of the fundamental degrees of freedom. In field theories, the variational approach is indispensable in 
discussing, for example, the underlying symmetries of the matter field, such as Noether's theorem, the quantization 
procedure, etc. Most of these formal aspects of the variational approach will not be much useful for hydrodynamical 
systems, except for the obvious symmetries required for the energy-momentum tensor. Thus one might find no 
point in discussing hydrodynamics from the action principle, once the equations of motion of hydrodynamics are well 
established in terms of the equations for the energy-momentum tensor. 

On the other hand, as is well-known, the variational approach has practical advantages besides its formal side. Once 
the variational principle is established, we can use the method to obtain the optimal parameters of a given family of 
trial solutions. 

In this paper, we derived the equations of motion of hydrodynamics starting from a very simple Lagrangian density. 
There it is seen that the roles of the continuity equation as a dynamical constraint and of the local thermodynamical 
relations are essential to arrive at the standard result of hydrodynamics. When the continuity equation is soluble, such 
a formulation in terms of the variational principle offers a powerful tool to obtain approximate solutions. For a system 
with a high degree of symmetry such as spherically symmetric system, we can establish the effective Lagrangian for the 
density profile function. Such an effective Lagrangian is quite useful for obtaining the approximate solutions for the 
hydrodynamical equation of motion in a simple manner. Even for the finite element discretisation of the hydrodynamic 
equation designed to a larger numerical solution, the variational approach may offer a physically optimized equation 
of motion, avoiding the mathematical instability with a relatively small number of degrees of freedom ||l^. As an 
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example of extreme simplified case, we apply the effective Lagrangian formulation for a gas bubble in a fluid and for 
the first time the relativistic version of the Rayleigh-Plesset equation is obtained. Such an approach will be useful 
for the analysis of the relativistic motions of blast waves in the models of gamma ray bursters Ijl^l , or the hot and 
dense droplet of QGP plasma, possibly formed in high-energy nuclear collisions. The application of our formalism for 
such processes is being planned. We have discussed also the introduction of pseudo-viscosity due to Neumann and 
Ritchmyer in the context of variational formulation. This will allow, for example, not only to treat relativistically the 
propagation of shock waves but also to introduce the finite relaxation time of turbulent flows in a phenomenological 
manner ||14| in the relativistic fluid dynamics. 

Our formalism will be useful in studying some problems of General Relativity, too. For a spherically symmetric 
system, a very simple Lagrangian density has been found. From this Lagrangian density we can show that all 
the known equations of the spherically symmetric system can be derived. We expect that, together with a good 
parametrization of the metric functions, approximate solutions to these otherwise difficult problems of stellar collapse 
or explosion with realistic equations of state can be obtained. Work on this line is in progress. 
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